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We study a classical model of Helium atom in which, in addition to the Coulomb forces, the 

radiation reaction forces are taken into account. This modification brings in the model a new 

qualitative feature of a global character. Indeed, as pointed out by Dirac, in any model of classical 

electrodynamics of point particles involving radiation reaction one has to eliminate, from the a priori 

conceivable solutions of the problem, those corresponding to the emission of an infinite amount of 

t^j- ' energy. We show that the Dirac prescription solves a problem of inconsistency plaguing all available 

ON ' 

£ — . models which neglect radiation reaction, namely, the fact that in all such models most initial data 

in 

lead to a spontaneous breakdown of the atom. A further modification is that the system thus 
acquires a peculiar form of dissipation. In particular, this makes attractive an invariant manifold 
of special physical interest, the zero-dipole manifold, that corresponds to motions in which no 

- i — i i 

energy is radiated away (in the dipole approximation). We finally study numerically the invariant 
measure naturally induced by the time-evolution on such a manifold, and this corresponds to 
studying the formation process of the atom. Indications are given that such a measure may be 
singular with respect to that of Lebesgue. 
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At variance with the existing works on the dynamical properties of classical 
atomic systems (see for example Refs. |lH8|), in this paper we take into account, 
for the case of the Helium atom, the effects on the dynamics due to the radiation 
emitted by the electrons during their motion. This introduces a qualitative 
change in the dynamics, with respect to the purely Coulomb model, because 
the model now corresponds to a non conservative dynamical system. This leads 
to the appearing of an "attractive" manifold of stable periodics orbits, where 
attractive is to be understood in a proper sense which will be explained below. 
The aim of this paper is to study the properties of such a manifold, through 
analytical and numerical methods. 



INTRODUCTION 



The studies on the dynamical properties of classical atomic models, in which only 
Coulomb forces are taken into account, have a long history, going back to Nicholson, 
Bohr and Langmuir (see Refs. lll-y). The main result found in such old works (see the 
work !ll| of Nicholson for the case of Beryllium, at those times called Nebulium) is that there 
exist periodic orbits having "normal mode frequencies" whose ratios are near to those of the 
observed spectral lines. Among the recent works on this subject a particularly relevant place 
is taken by the papers of De Luca jl, 8] on the Helium atom. This author went beyond the 
Coulomb approximation, taking retardation of the electromagnetic forces into account. This 
was obtained through a perturbation scheme which, truncated in a suitable way, leads to 
the conservative Darwin lagrangian that one finds illustrated in common textbooks 
In such a way the scale invariance of the purely Coulomb model was eliminated, and the 
result of Nicholson could be improved, showing that for the Helium atom there exist periodic 
orbits leading to dynamical frequencies which themselves, and not only their ratios, have a 
rather good agreement with the empirical spectral frequencies. 

On the other hand, it was known since the time of Nicholson, and confirmed for example 

fl 

by Poirier 14], that all such periodic orbits are unstable. More recently Yamamoto and 
Kanekojsl, IJ showed, for the case of Helium, that the instability of the periodic orbits 
actually corresponds to a much more general and acute form of instability. Indeed they 
found that the vast majority of initial data with negative energy lead to the autoionization 
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of the atom, namely, to motions in which one of the electrons is expelled. One can see that 
this occurs also for a generic atom, and even in a stronger way, inasmuch as one meets with 
motions in which all electrons but one are expelled. So, the whole classical theory of atomic 
models seems to be plagued by a general failure, for which no remedy is known. 

One usually takes the pragmatic attitude of just forgetting the autoionizing motions. One 
should however indicate an internal dynamical justification of a more general character for 
such a selection of the initial data. Moreover, one might hope that such a general prescription 
would also automatically explain why the system chooses to fall (through some peculiar kind 
of dissipation) on the physically interesting periodic orbits, for which an agreement betweeen 
dynamical and spectral frequencies is found. Obviously, satisfying the latter point requires 
to abandon a purely lagrangian, and thus conservative, approach to the problem, by taking 
electromagnetic radiation into account. 

This is precisely what we do in the present paper, where we take into account the emission 
of radiation due to the accelerations of the electrons in the simplest possible way, namely, 
by adding to the Coulomb forces the radiation reaction ones. So, in a sense our model is 
complementary to that of De Luca, who introduces only the conservative terms produced 
by the Darwin approximation scheme, and neglects dissipation. We leave for future work 
the study of a model which takes both features into account. 

It will be shown here that the modification of the Coulomb model that takes radiation 
into account, in the first place gives a model in which the autoionization problem no more 
arises. This is due to the fact, first pointed out by Dirac in the general context of classical 
electrodynamics of point particles 12], that most initial data in the phase space suited to the 
model lead to motions in which an infinite amount of energy is emitted, so that the definition 
of the model has to be complemented by the explicit prescription that such initial data have 
to be discarded. The remaining initial data (constituting a set that will be called here the 
Dirac or physical manifold) will be shown to lead to motions in which the phenomenon of 
autoionization no more occurs. So, the elimination of the initial points in phase space leading 
to autoionizing solutions appears no more as a special trick to be strangely introduced ad 
hoc, but rather as a particular case of a completely general prescription that, following 
Dirac, has always to be introduced when dealing with classical models of matter-radiation 
interaction involving point particles. 

Then we study the dynamics on the Dirac physical manifold, in which the Dirac pre- 



cription turns out to introduce a dissipation af a peculiar type. We show that the Dirac 
prscription allows one to find (and in a easy way) all periodic orbits, proving furthermore 
that they form in phase space a manifold. This is just the zero-dipole manifold, which is 
constituted of phase-space points leading to motions that do not radiate energy away in the 
dipole approximation. A major aim of this paper is to study such a manifold. We first show 
that it is an attractor, for initial states with negative energy. Then we study, by numerical 
methods, the invariant measure naturally induced on it by the time-flow. This amounts to 
studying the formation process of the Helium atom, namely, motions having initial data with 
the two electrons coming from infinity that are then captured by the nucleus, and finally 
fall on the zero-dipole manifold, where emission of energy comes to an end. 

From the point of view of the theory of dynamical systems, an interesting result we find 
is that, while the attractor (the zero-dipole manifold) is really simple as a manifold (being 
just a portion of a linear subspace in the system's phase space), it is the invariant measure 
induced on it by the time evolution that is very peculiar. Indeed it is presumably singular 
with respect to the restriction of the Lebesgue measure, possibly having a fractal structure. 

In Section 2 the model is introduced, and a preliminary analytical discussion of its prop- 
erties is given. In Section 3 the numerical method for obtaining the invariant measure is 
described, and the numerical results are presented. The conclusions follow. 



II. THE MODEL 

A. The radiation reaction for a single point charge 

Let us recall that, in the case of a single point-charge, the emission of radiation is taken 
into account without introducing the infinitely many degrees of freedom of the electromag- 
netic field, through the expedient of adding in the equation of motion for the particle an 
"effective force". This force, which is traditionally called radiation reaction force and de- 
noted by K. rad , is given by 

K^=jgx, (1) 

where e is the charge of the particle, c the speed of light, and x the position vector of 
the particle. The procedure which leads to such a force goes back to Planck, Lorentz and 



Abraham, and was given a final form in the work of Dirac 



121 ] of the year 1938, in which 



an extension to the relativistic case was performed in an extremely elegant way. In its most 
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elementary form, which is illustrated in standard textbooks (see Refs 
the procedure amounts to requiring that the effective force K md produces an energy loss 
consistent with the power radiated away by an accelerating particle according to the Larmor 
formula, namely, (2/3) (e 2 /c 3 ) |x| 2 . As is well known, this force can also be interpreted as 
being produced by the regular part of the self-field of the particle, the divergent part having 
been reabsorbed through mass renormalization. 

We illustrate now, in the simple case of the free particle in the nonrelativistic approxi- 
mation, how the "runaway" solutions then show up. In terms of the particle's acceleration 
a = x, the equation of motion for the free particle takes the form 



with 

2 e 2 



(2) 



3 mc 3 

(m being the particle's mass), with solutions a(£) = a exp(t/e) depending parametrically on 
the initial acceleration a . So, for generic initial data the solution exponentially explodes, 
i.e., presents runaway character, a fact which makes no sense for a free particle. Thus 



Dirac [12|] introduces the explicit prescription that, for the case of a free particle, only the 
solutions with vanishing initial acceleration, ao = 0, be retained. Now, the phase space 
V mathematically suited to the considered third-order equation of motion, is the vector 
space R 9 , a point of which is defined by the coordinates (x, v, a) of position, velocity and 
acceleration. In such an ambient phase space V, the physically meaningful phase space is 
thus the subset characterized by motions of nonrunaway type. This we will call "Dirac or 
physical manifold" , and coincides with the hyperplane a = 0. 

More in general, in the presence of an external force vanishing at infinity, one should 
assume that for scattering states, in which the particle eventually behaves as a free one, the 
physical or nonrunaway manifold be defined by the asymptotic condition 

x(t) ^ for t +oo . (3) 

Prescriptions of such a type may be formulated in a different way, which captures another 
side of the problem. One makes reference to the total energy emitted by the particle, namely, 
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the quantity 

AE = me / |x(t)| 2 dt , (4) 
Jo 

and a physically natural requirement is then to restrict oneself to motions for which the 
amount of energy radiated away during the whole motion is finite, i.e., one has AE < oo. 
For smooth motions, such as those we are considering which are solutions of an ordinary 
differential equation, the latter condition implies the asymptotic condition ([3]). Thus con- 
dition ([3]) turns out to be justified also for motions not having a scattering character (as 
the bound states), for which the argument referring to the eventual free-type motion of the 
particle does not apply. 

The asymptotic condition ([3]) will be referred to as the Dirac prescription. Due to the 
global character (with respect to time) of the latter, the dynamics on the Dirac physical 
manifold thus acquires quite peculiar fetures [16]; for example, it can be folded 17, 18]. 



We will give arguments indicating that, for the Helium atom model with radiation re- 
action, the analog of the Dirac prescription automatically eliminates the autoionization 
problem, and thus overcomes the main inconsistency of all classical models which neglect 
radiation. 



B. Definition of the model 

In the case of several particles one can operate as for just one particle, by introducing 
a suitable radiation reaction force acting on each particle, in such a way that the power 
dissipated by the whole system be equal to that emitted in the dipole approximation (see 
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sec. 9.2). In the case of two particles the radiated power is me(xi + x 2 ) 2 , where x 1; x 2 
are the position vectors of the particles (here, the position vectors of the two electrons with 
respect to the nucleus, assumed to be fixed at the origin of an inertial frame). This gives 
for the effective forces K™ d and K™ d acting on each electron the same expression, namely, 
mefxi +x 2 ). Thus the model is described by the system of equations 

2e 2 xi e 2 xi - x 2 ... . 

xi = H 3 h £ (xi + x 2 ) 

m 77 m r^2 

2e 2 x 2 e 2 x 2 - xi ... 

x 2 = 3-H 3 h£(x 2 + xi) (5) 

m r 2 m r\ 2 

(in the Gauss system, with e given by (j2J)). 



One obtains the same expression for the radiative force acting on each particle, also if 
one considers it as due to the sum of the (regular part of the) self-field, and of the retarded 
field produced by the other charge, expanded with respect to the distance. One can check 
this fact through computations which are completely standard, though a little cumbersome 
and not particularly illuminating. So we omit them. The only point worth of mention is 
that the expansion in the distance should be performed up to third order and not just to 
the second one as one finds in the textbooks; this is indeed the point which is responsible 
for the appearing of the relevant nonconservative term proportional to the third derivative 
of the center of mass. 

System fl5]), when expressed in terms of the center of mass x cm and of the relative position 
vector r, defined by 

x cm = (xi + x 2 )/2 , r = x 2 - xi , 



takes the form 



2e r 2e / x cm r/2 x cm -|- r/2 



s + 



m r° m 



3 yyi \ I I r\ I ^ / c\ I ^ 



^cm ^ / 2 P^cm ^/2 



o •■- •■ \ ^ f Xcm r/2 x cm H~ r/2 \ 

zeyicm - x cm + — y- - 3 + - ,3 ) » W 

in which there appears only one third-derivative term, x cm , which enters just the equation 
for the center of mass. In conclusion, when radiation reaction is taken into account, the 
equations of motion in the variables x cm , r are exactly the same as for the purely Coulomb 
model, with the only addition of the radiation reaction force in the equation for the center 
of mass. 

So the ambient phase space V suited for the motions of the two point electrons is the 
vector space R 15 , a point of which is defined by the coordinates 



Now, the radiation reaction term acts in the equation for the center of mass just as in the 
case of a single particle. Thus, in the spirit of the classical work of Dirac, we will define our 
model by complementing the system of equations (J5j) or ([6]) through the Dirac prescription: 
in the ambient phase space V = R 15 the only admitted points are those leading to motions 
that satisfy the asymptotic condition 

x cm (t) ->■ for f-^ +00 . (7) 



As in the case of the single particle, this corresponds to the physical condition that 
the amount of energy radiated away during the whole motion be finite. This prescription 
implicitly introduces a selection among the allowed initial data, and the subset of the allowed 
points in the ambient phase space V will be called the Dirac or physical manifold. We will 
show that this restriction eliminates the autoionization problem, and makes the zero-dipole 
manifold become attractive for initial data of negative energy. 



C. Simple analytical deductions 

One immediately sees that there exists a particular invariant submanifold of the Dirac 
physical manifold, which will play a fundamental role in this work. It is the zero-dipole 
manifold, defined by the condition x cm (t) = for all times, which corresponds to the hyper- 
plane 

^cm , x cm , x cm . (8) 

Invariance is immediately checked. One also easily checks that such a manifold is composed 
of orbits which are solutions of the purely Coulomb model in the unknown r(t) with a 
suitable charge, and are thus periodic in the case of negative energy. 

We will make use of an energy theorem. This is immediately established by multiplying, 
as usual, equations §5§ by x x and x 2 respectively, adding them and using the Leibniz formula 
for the product 'x cm x cm . This gives 

where T is the kinetic energy, V the potential energy corresponding to the Coulomb forces, 
and 4me x cm ■ x cm the so-called Schott term. So, the quantity 

£ = T + V - Amex crn ■ x cm 

may be simply called the energy of the system, (while the quantity E = T + V may be called 
the mechanical energy), and £ turns out to be a non increasing function of time. Obviously 
the two energies coincide, £ = E, on the zero-dipole manifold. 

From the energy theorem one immediately gets that periodic orbits necessarily belong 
to the zero-dipole manifold. Indeed, just integrating the energy equation over a period, a 
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periodic orbit is seen to necessarily have vanishing center of mass acceleration, and con- 
sequently is seen to belong to the zero-dipole manifold. In addition, the portion of the 
zero-dipole manifold with negative energy is wholly foliated by periodic orbits which, as 
already mentioned, are just the periodic orbits for the two-body Coulomb problem with a 
suitable charge. 

Using the energy theorem one can also make more precise the notion of an autoionizing 
motion. We recall that the instability problem that plagues the purely Coulomb model is 
that the vast majority of initial states with negative energy lead to motions in which one of 
the electrons is expelled to infinity, so that the atom is unstable. We show instead that the 
analogous situation does not occur in the model with radiation reaction, for initial states 
with negative energy. 

Indeed, if autoionization occurs with the two electrons both escaping to infinity, then 
both the potential and the Schott term finally vanish, so that the energy finally becomes 
positive, i.e. it has increased, against the energy theorem. If instead autoionization occurs 
with only one electron, say the first one, escaping to infinity (in a nonrunaway fashion), then 
the equation of motion for x 2 reduces to 

2e 2 x 2 

x 2 = T + ex 2 , (9) 

m r'2 

which is the equation for just one electron with radiation reaction, in the external field of 

the nucleus. On the other hand, as shown in papers Q and 13, starting from initial data 

with negative energy, equation fl9]) admits only runaway solutions. 

Furthermore, the zero-dipole manifold is an attractor for initial states with negative 

energy, as a consequence of the asymptotic Dirac condition ([7j). This can be seen by rewriting 

the equation for the center of mass in the integro-differential form 

^ r+00 

*cm{t) = --p^ £ J e~ s/2£ f(x cm {s),r{s))ds 

where / is the function defined by the last two terms at the right hand side of the second 
equation in (JBJ). So, from the de L'Hopital rule the property / — >■ for t — > +00 follows. 
Finally, the result x cm (t) — > for t — > +00 follows from the form of the function /, using 
the previously proven fact that for negative energies the motion of the center of mass is 
bounded. 

In conclusion, following the Dirac prescriprion of restricting the ambient phase space to 
the Dirac physical manifold, on the one hand the difficulty of the generic autoionization is 



eliminated, and on the other hand our dynamical system presents a dissipative character, 
because all initial states with negative energy are definitively attracted to the invariant 
zero-dipole manifold, corresponding to motions that do not radiate energy away (in the 
dipole approximation). However, in the next section we will show that the electrodynamical 
dissipation due to radiation reaction has a peculiar character with respect to that of the 
familiar dissipative systems such as the Lorenz one. 



III. THE NATURAL INVARIANT MEASURE ON THE ATTRACTOR: NUMER- 
ICAL RESULTS 



It is simple to check that the Lebesgue measure in phase space, restricted to the zero- 
dipole manifold, is invariant under the flow. However, other, actually infinitely many, in- 
variant measures exist, this being a general property of any continuous group $* of diffeo- 



morfisms on a smooth manifold Ai. Indeed, following Krylov and Bogolyubov [21], for any 
given measure /io with fi (Ai) < oo, a corresponding invariant measure \i can be constructed 
by defining 

1 /"* 

fj.(B)= lim - / // ($~ s £)ds. (10) 

t-++oo t J 

The intuitive meaning of the measure /i is the fraction of the initial points of phase space 
(chosen according the given probability /Jo) which fall on any given set B. So, if we choose 
for /io the restriction of the Lebesgue measure to the Dirac physical manifold, and take initial 
data with the two electrons coming freely from infinity, then from the physical point of view 
the measure /i will describe the state of the atom at the end of the formation process. 

Our aim is now to find a numerical scheme for constructing the measure fi. For the sake 
of simplicity, in the implementation we will restrict ourselves to the case of planar motions. 
So, the phase space corresponding to equations ([6]) has dimension 10, and the Dirac physical 
manifold has at most dimension 8, while the zero-dipole manifold has dimension 4. 

It goes without saying that the motions are to be obtained by integrating the system 
backward in time. Indeed the Dirac physical manifold is defined only implicitly, through the 
asymptotic condition ([7]). Thus one concretely has to operate in the ambient phase space V, 
in which the Dirac manifold has vanishing measure. Besides, even if one were able to take 
initial data on the Dirac manifold, by integrating forward in time numerical errors would 
make the orbits escape from it, and in a runaway fashion. Such problems are overcome by 
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the standard procedure of integrating the equations of motion backward in time. In such 
a way one is guaranteed to approach (in a time of the order 2e) the Dirac manifold. By 
the way, the integration beackward in time is just the one needed in order to compute the 
asymptotic measure fi according to the definition ffTUl) . 

From the numerical point of view, the only practical way in which a measure can be 
described, is to compute its density p(x). Here we denote by x the coordinates of a point 
in the ambient phase space V. So we implemented the following algorithm. Taking x on 
the zero-dipole manifold as initial datum, we follow, backward in time, the time evolution 
of a small hypercube of equal sides \Sxj\ (j = 1, . . . , 10), having x as one of its vertices, 
i.e., we follow the evolution of the points Xj = x + 5x.j, and define, as usual, the side of 
the evolved hypercube by 8xj(t) = — $*x. Notice that the hypercube does not belong 
to the zero-dipole manifold (not even to the Dirac one); however, in our specific case, the 
integration of the equations of motion, backward in time, insures that the trajectories fall 
on the Dirac manifold after a very short transient (of the order 2e), because the runaway 
component of the motion is exponentially damped as one proceeds backward in time. So, the 
hypercube becomes squeezed on the physical manifold, and one of its 8-dimensional faces 
becomes tangent to the manifold. To get the density p(x), it is then sufficient to compute 
the area of the latter face, divided by the initial area, and then calculate its time-average 
along the trajectory of x (which obviously belongs to the zero-dipole manifold for all times). 

We chose to numerically integrate the system with a fourth order Runge-Kutta method 
and autoadaptation of the time grid. We take units with m — 1 and e = 1, completing 
them by taking as unit of length the Bohr radius Rb, as usually done in atomic physics. 
In this way e would take the value 1 / (137) 3 . However, as we are only interested in the 
qualitative behaviour of the system, in order to reduce the computational time we decided 
to take instead e = 10~ 2 . The autoadaptation of the time grid was implemented by choosing 
the integration step h according to 



where ho = 10~ 4 , while r\ , r-i and r\2 are the distances of the first and second electron from 
the nucleus and their mutual distance respectively, while r min = 0.001 and r max = 1 are 
factors introduced to avoid that the integration step becomes too small or too large. The 
sides \5xj\ of the initial hypercubes were all taken equal to 10~ 4 . 



h = ho ■ max 




11 



Once the trajectories are known (numerically), there remains the problem of how to 
identify the face of the hypercube lying on the Dirac manifold (which is utterly unknown). 
In order to solve this problem, one can start from the fact that all sides Syij are essentially 
contained in the tangent plane to the manifold. Then this plane can be found by determining 
the 8-dimensional hyperplane to which the vectors 5xj are closer. As an hyperplane is 
defined by two independent unit vectors orthogonal to it, say ai and a 2 , then the tangent 
plane will be defined by the unit vectors which minimize £V (a x • 5x.j) 2 + (a 2 • 5xj) 2 , i. 
e., the sum of the distances of the points 5x.j from the plane. In other terms one has to 
minimize the function 

S(a ± , a 2 ) = f ajBsix + a 2 r 5a 2 
where the symmetric matrix B is defined by 

B = [5xi, . . . , <5xi ] T • [5xi, . . . , 5xio] . 

As is well known, the two independent unit vectors which minimize S are the eigenvectors 
corresponding to the two smallest eigenvalues Ai, A 2 of B. With a little reflection, one 
also gets that the area of the tangent face is nothing but the square root VA3 • • • A 10 of the 
product of the remaining 8 eigenvalues Aj, i — 3, . . . , 10 of the matrix B. 

In summary, the numerical procedure consists in taking for p(x) the time-average of 



VA3 . . . Aio- Notice that, as occurs in the computation of the Lyapunov exponents, the sides 
Sx.j(t) can grow too much during the evolution. If, for some j, \5x.j\ becomes larger than 
10~ 3 , we renormalize the cube choosing a new one with sides directed as the eigenvectors of 
B and again of size 10~ 4 , and the average is then accordingly modified. 

Just to give an indication of the form of the natural measure, in the zero-dipole manifold 
we chose 10 5 initial points x uniformly distributed in the cube 

r G [0.5,4.0] x [0.5,4.0] 
r G [-1.5,1.5] x [-1.5,1.5] , 

with the condition that they have negative energy. We checked that averaging over a period 
is enough to insure that the time average is stabilized. 

Now, the zero-dipole manifold is 4-dimensional, so that in principle p depends on four 
variables (i.e., position and velocity). Recall however that on such a manifold the motion 
is a purely Coulomb one, so that it is integrable. Then, thinking in terms of action-angle 
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FIG. 1. Invariant density p on the invariant zero-dipole manifold, projected on a rectangle of the 
(E, M) plane, of energy and angular momentum. 
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E 

FIG. 2. Same as Figure [TJ with the values of p represented in a gray decimal logarithmic scale, 
covering a thousand orders of magnitude. 




FIG. 3. Same as Figure [H restricted to a smaller portion of the original rectangle, with the 
ordinates in a decimal logarithmic scale as in Figure [2] . 
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variables, the interesting ones are just the two actions, which in our case are E' 1 / 2 , E 
being the (mechanical) energy, and the angular momentum M, whereas the angles, i.e., the 
orientation of the orbit and the starting point of the motion on the orbit, are irrelevant. So we 
chose to integrate away the two angles from the density p, and to analyze the corresponding 
reduced density, p(E, M) (by abuse of language still denoting it by p). 

The result is shown in Figures [TH3J In Figure [1] the reduced density p is reported versus 
E and M. While a priori the measure could be spread over the whole rectangle, the figure 
shows that it is actually concentrated in an extremely small area. This is more clearly 
exhibited in Figure [2J which still covers the whole rectangle, but gives the density in a 
decimal lograthmic gray scale, which covers a thousand orders of magnitude. Figure [3] is 
instead the same as Figure [TJ restricted to a smaller portion of the rectangle, still with 
the ordinates in a decimal logarithmic scale. It shows that the function p(E, M) has a very 
complex structure, with values ranging over a thousand orders of magnitude in an apparently 
nonsmooth way. This fact might suggest that the measure is not absolutely continuous with 
respect to the restriction of the Lebesgue one, possibly having some kind of fractal structure, 
but we leave this problem for possible future work. 

In any case, we have shown that the peculiar character of dissipation entering electrody- 
namics of point particles restricted to the Dirac physical manifold, makes the final dynamics 
quite different from that of more familiar dissipative systems such as the Lorenz one. Indeed 
in the latter case it is the invariant attractor that has a strange character, while the natural 
invariant measure induced on by the dynamics is somehow trivial, as corresponding to that 
of a hyperbolic system. In the former case, instead, the invariant manifold is trivial (being 
just the hyperplane x cm = 0,x cm = 0,x cm = in the ambient phase space), while it is the 
invariant measure which apparently has a strange character. 

IV. CONCLUSIONS 

We investigated the main modifications that are introduced in the Coulomb Helium atom 
model when the radiation reaction forces are taken into account in the simplest possible way, 
and the asymptotic Dirac prescription is consequently introduced, according to which the 
phase space has to be restricted to the submanifold of points leading to motions with a finite 
emitted energy. 
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The first qualitative result we found is that the problem of the breakdown of the atom 
(corresponding to autoionization for generic initial data with negative energy), which makes 
the Coulomb model inconsistent, is now eliminated. The second one is the existence of an 
invariant manifold of stable periodic orbits, that moreover is attractive. Such a manifold 
is the zero-dipole one, on which the system does not radiate energy away, in the dipole 
approximation. 

Finally, the invariant measure naturally induced by the time-flow on the zero-dipole 
manifold was studied numerically. We showed that such a measure is far from trivial, 
being presumably non absolutely continuous with respect to the restriction of the Lebesgue 
measure, possibly with some kind of fractal structure. This is at variance with the known 
examples of dissipative systems, for which the measure is trivial, while it is the attractor 
that is strange (i.e. a fractal set). This is apparently a consequence of the peculiar character 
of the radiation reaction force, which makes the system dissipative when restricted to the 
Dirac physical manifold, while making it expansive in the ambient phase space. 

A very interesting problem that remains open is that of comparing the special periodic 
orbits which are selected on the zero-dipole manifold as the most probable ones according 
to the natural measure, with those, apparently of special physical relevance, studied by 
De Luca. We leave this point for future work. 
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